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Abstract It is described how one comes to the Wigner-Dyson random matrix the- 
ory (RMT) starting from a model of a disordered metal. The lectures 
start with a historical introduction where basic ideas of the RMT and 
theory of disordered metals are reviewed. This part is followed by an 
introduction into supermathematics (mathematics operating with both 
commuting and anticommuting variables). The main ideas of the su- 
persymmetry method are given and basic formulae are derived. Both 
level-level correlations and fluctuations of amplitudes of wave functions 
are discussed. It is shown how one can both obtain known formulae 
of the RMT and go beyond. In the last part some recent progress in 
the further development of the method and possible perspectives are 
discussed. 
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Introduction 

Wigner-Dyson Theory 

According to basic principles of quantum mechanics the energy spec- 
trum of a particle in a limited volume is discrete. The precise values of 
the energy depend on the boundary conditions and the interactions in 
system. In many cases these quantities can be calculated with a certain 
accuracy. However, often the interactions are so complicated that cal- 
culations for the single levels become impossible. On the other hand, 
the complexity of the interactions can lead to the idea of a statistical 
description in which information about separate levels is neglected and 
only averaged quantities are studied. Density of states, energy level and 
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wave functions correlations are quantities that can be studied in the 
statistical approach. Sometimes, it is sufficient to study the average of, 
e.g., the density of states or its variance. In other cases one may be in- 
terested in a full statistical description that can be achieved calculating 
distribution functions. Studying the level statistics is, to some extent, 
analogous to the statistical study of the motion of atoms and molecules, 
which is the subject of statistical physics. 

The idea of the statistical description of the energy levels was first 
proposed by Wigner [1] for study of highly excited nuclear levels in 
complex nuclei. In such nuclei a large number of particles interact in 
an unknown way and the main assumption was that the interactions 
were equally probable. Of course, in order to specify the meaning of the 
words "equally probable" one had to formulate a statistical hypothesis 
in terms of a probability distribution that would play the role of the 
Gibbs distribution. 

This was done in Ref. [1] in the following way. Choosing a complete 
set of eigenfunctions as a basis, one represents the Hamiltonian H as a 
matrix with matrix elements H mn . The matrix elements H mn are as- 
sumed random with a certain probability distribution. It is clear that 
the distribution should not depend of the basis chosen, which implies an 
invariant form of the distribution function. In the language of the ran- 
dom matrices H mn the corresponding distribution function can contain 
only Trf (H), where / is a function. 

The first statistical theory [1] was based on a Gaussian distribution. 
According to the Gaussian statistical hypothesis a physical system hav- 
ing N quantum states has the statistical weight D (H) 



In Eq. (1) the parameter a is a cutoff excluding strong interaction, H is 
a random N x N matrix, and A is a normalization coefficient. 

It is important to emphasize that the weight D (H) is rather arbitrary 
and other forms for the distribution functions can be suggested. Of 
course, the dependence of the mean energy level spacing A (e) 



on the energy e is different for different distributions D (H). (In Eq. (2) 
the symbol stands for the averaging with the distribution D (H)). 

For example, the distribution function D (H) for the Gaussian distribu- 
tion, Eq. (1), has the form of a semicircle (Wigner semicircle law) (for 
a review, see, e.g. Ref. [2]). 




m,n=l 
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What is more interesting, level correlations described, for example, by 
the level- level correlation function R (ui) 



prove to be universal in the limit N — > oo provided the energy uj is much 
smaller than the characteristic scale of the variation of A (e) . For the 
Wigner semicircle law, the latter condition means that the energy e is 
not close to the points eo, — £q, where the quantity A -1 (e) proportional 
to the average density of states turns to zero. 

In the absence of magnetic interactions violating the time reversal 
symmetry, the wave functions and matrix elements H mn in Eq. (1) can 
be chosen real. In this case the statistical properties of the systems are 
described by real symmetric random matrices. The ensemble of real 
symmetric matrices with the Gaussian distribution is often called the 
Gaussian orthogonal ensemble (GOE). 

If the magnetic interactions are present in the system, the time rever- 
sal and spin-rotation symmetry is violated and the wave functions are 
no longer real. This means that one should deal with general Hermitian 
matrices without any additional symmetry and integrate over the matrix 
elements H mn using only the constraint H mn = (H nm )*. This system is 
called the Gaussian unitary ensemble (GUE). 

The third possible type of the symmetry arises when the system is 
time-reversal invariant but does not have central symmetry. In this case 
it is also impossible to make all the matrix elements real. Nevertheless 
an additional symmetry exists in this case. According to the Kramers 
theorem all levels of the system remain doubly degenerate and every 
eigenvalue of the matrix H must appear twice. Matrices consisting of 
real quaternions H mn of the form 



and satisfying the condition H mn = (H nm ) have this property. The cor- 
responding ensemble is called the Gaussian symplectic ensemble (GSE). 

Somewhat different distribution functions were introduced later by 
Dyson [3] who suggested characterizing the system not by its Hamil- 
tonian but by an unitary N x N matrix S whose elements give the 
transition probabilities between states, where, again, N is the number 
of the levels. This matrix is related to the Hamiltonian H of the system 
in a complicated way that is not specified in the theory. According to 
the Dyson hypothesis the correlation properties of n successive energy 
levels of the system (n <C N) are statistically equivalent to those of n 



R (u) = A 2 (e) (Tr5 (e - H) Tr5 (e - u - H)) 



(3) 
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successive angles provided all the unitary matrices S have equal proba- 
bilities. Again, depending on the symmetry, one can distinguish among 
the three different ensembles. The corresponding ensembles are called 
Circular ensembles. 

As concerns the complex nuclei, the orthogonal ensembles are most 
relevant for their description because in order to change the level statis- 
tics one needs, e.g. huge magnetic fields that hardly exist. However, 
the other two ensembles have been under intensive discussions for prob- 
lems of mesoscopic physics (for a review, see [4]). Moreover, it has been 
realized that one might formulate additional ensembles of e.g. chiral 
matrices relevant for studying properties of models for QCD (for a re- 
view, see the lecture by Jac Verbaarschot [5] ) . It has been proven later 
that in total 10 different symmetry classes exist [6]. Most of the new 
ensembles may be relevant to different disordered mesoscopic systems 
due the presence of, e.g., superconductivity or additional symmetries of 
the lattice. 

I do not plan to discuss in these lectures the non-standard symmetry 
classes and restrict myself by the Wigner-Dyson (WD) statistics. It is 
relevant to say that calculation of the level-level correlation function, Eq. 
(3), starting from the Gaussian or Circular ensembles is not a simple 
task. The conventional method of the evaluation is using orthogonal 
polynomials [2]. The procedure is not difficult for the unitary ensembles 
but one has to put a considerable effort to perform the calculations for 
the orthogonal and symplectic ones. As it has been mentioned, the final 
results are universal in the limit N — > oo and can be written for the 
orthogonal, unitary and symplectic ensembles in the form 

, . sin 2 £ d f smx\ f°° sinxt , , A . 

iWM = i- — — / — * (4) 



D f \ -t S[li2x fti\ 
Runit M = 1 — (5) 

sin 2 x d ( sinx\ f 1 sinxt 
R sympl (w ) = l-^ + -(_jjf —eft (6) 

where x = ttlu/A, and A is the mean level spacing. 

The functions R(co), Eqs. (4, 5, 6), tend to 1 in the limit u> — ► oo, 
which means that the correlations are lost in this limit. In the opposite 
limit x — > they turn to zero as x 13 , where (3 = 1,2,3 for the orthogonal, 
unitary and symplectic ensembles, respectively. This means that the 
probability of finding a level at the distance oj from another level decays 
at small to. The effect is know as "level repulsion". It is important that 
Eqs. (4-6) describe the level-level correlation function, Eq. (3), for both 
the Gaussian and Circular ensembles. 




Figure 1. Structure of a granular metal 



Small disordered particles 

Nowadays the relevance of the random matrix theory (RMT) to meso- 
scopic physics is almost evident and it is the starting point of many 
works on transport in quantum dots, electromagnetic response of metal- 
lic grains, etc. However, it took quite a long time before the ideas of the 
RMT penetrated from nuclear to condensed matter physics. 

There were several reasons for this slow development. First, until the 
end of 60's of the last century most of the objects studied in condensed 
matter physics were macroscopic and the discreteness of the energy lev- 
els could be neglected. Therefore the question about the level statistics 
was not so interesting. Second, from the theoretical point of view it was 
not clear at all how one could come to the Wigner-Dyson level statistics 
starting from the Schrodinger equation. It was clear that making per- 
turbation theory in interaction or disorder did not lead to the anything 
that would resemble Eqs. (4-6). 

As soon as experimentalists started investigation of granular materi- 
als (which happened in 60's), a theory that would describe small metal 
systems became nevessary. Such materials consist of small metal parti- 
cles (grains) with the diameter down to 10 — lOOA. These grains can be 
covered by an insulator and therefore be well isolated from each other. 
A schematic view of the pattern can be found in Fig.l 

It is clear that one can speak now about discrete levels and study 
their statistics. Of course, one needs low temperatures in order to pre- 
vent inelastic processes smearing the levels but this is not so difficult 
(temperatures < IK can be sufficient). 

In practice, the form of the grains can be not very regular, they 
can contain defects and impurities and therefore the energy spectrum 
strongly fluctuates from grain to grain. All this true even if one neglects 
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the electron-electron interaction. So, one naturally comes to the idea to 
describe the levels statistically. 

The first work on the application of the RMT to small metallic grains 
was done by Gorkov and Eliashberg (GE) [7]. These authors studied 
the electromagnetic response of the system of the grains and therefore 
they needed an information about the level-level correlation in a single 
grain. As the RMT is purely phenomenological and nothing is assumed 
about the origin of the randomness, this theory was taken by GE to 
describe the correlations. Starting from the explicit form of the level- 
level correlation function R(uj), Eqs. (4-6), they calculated the desired 
physical quantity. 

GE identified correctly physical situations when the three symmetry 
classes might be used. In the absence of any magnetic and spin orbit in- 
teractions they suggested to use the orthogonal ensemble. If a magnetic 
field is applied or there are magnetic impurities in the grains, the unitary 
ensemble should be applicable. If there are no magnetic interactions but 
spin-orbital impurities are present, the grains should be described by the 
symplectic ensemble. 

Being the first application of the RMT in solid state physics, the paper 
[7] remained the only application during the next 17 years. This is not 
surprising because using calculational schemes existed in that time no 
indication in non-trivial level correlations could be seen. Let us discuss 
this point in more details. 

Studying a disordered system and neglecting electron-electron inter- 
actions one can start with the following Hamiltonian 

H = e(p) + U(r) (7) 

where e (p) is the operator of the kinetic energy and calculate the Green 
function G £ = (e — H)~ l performing an expansion in the disorder po- 
tential U(r). Usually, it is assumed that the U (r) is random and its 
fluctuations are Gaussian with 

<£/(r)) =0, (U(r)U(r>)) = ^-6(r-r>) (8) 

where v is the density of states and r is the mean scattering time. 

Making the perturbation theory in the random potential U (r) and 
averaging over this potential can be done using the "cross technique" [8]. 
Only diagrams without intersections of the impurity lines are important 
in the standard approximation of the weak disorder. For the one particle 
Green function the typical diagrams are represented in Fig. 2. 
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Figure 2. Typical diagrams for the Green function in a disordered metal 

As a result, one obtains for the averaged Green function (G) the 
following expression 

° = e - e (p) ± i/2r (9) 

where "+" corresponds to the retarded and "— " to the advanced Green 
functions. 

Standard calculations based on a summation of "ladder diagrams" 
lead in this case to the classical Drude formula for, e.g., the conductivity 
a(u) 

a {oj) = a ° , a = 2e 2 vD Q (10) 

1 — IUJT 

where Dq = Vqt/3 is the classical diffusion coefficient, uj is the frequency 
and e is the electron charge. 

So long as the grain size remained larger than the atomic distances, no 
deviations from this formula could be found and therefore the question 
about the applicability of the RMT remained open. It is relevant to 
mention works by Kubo performed in approximately the same time. 
In Ref. [9] he argued that even very small irregularities (with size of 
the order of atomic distances) of the shape of the metallic grains must 
lead to lifting of all degeneracies of eigenstates that are present in ideally 
spherical particles. This lead him to the conclusion that the mean energy 
level spacing A is inversely proportional to the volume V of the particle 

A = (vV)- 1 (11) 

where v is density of states at the Fermi surface of the metal. Later Kubo 
suggested [10] that the spacing distribution had to follow the Poisson 
law. The latter differs essentially from the Wigner-Dyson RMT by the 



Figure 3. Cooperon: a singular correction to conductivity 



absence of the level repulsion. The second work by Kubo, Ref.[10], was 
published several years after the work by Gorkov and Eliashberg, Ref. 
[7], and this shows that the applicability of the Wigner-Dyson statistics 
to the small metal particles was far from being established in that time. 

The situation started to change only at the end of 70's with the new 
developments in the theory of Anderson localization. In the publication 
[11] a new scaling idea was put forward for description of disordered sam- 
ples of an arbitrary dimensionality. The most unusual was a prediction 
that two dimensional initially metallic samples could not remain metals 
in the presence of an arbitrary weak disorder and had to acquire insu- 
lating properties. Again, using Eq. (10) it was not clear why something 
had to happen in two dimension (the localization in one-dimensional 
chains had been proven before and it was clear that the ladder diagrams 
summation leading to Eq. (10) was not sufficient for that case). 

Trying to understand how something unusual could happen in two di- 
mensions Gorkov, Larkin and Khmelnitskii [12] investigated more com- 
plicated diagrams and found that a certain class of diagrams could lead 
to a divergence in any dimensionality d ^ 2. These are "fan" diagrams 
with the maximal number of crossings. They are represented in Fig. 3 
and their sum is a new effective mode that is usually called "cooperon" . 
This mode has a form of a diffusion propagator and its contribution to 
the conductivity can be written in the form 



It is clear from Eq. (12) that in the limit of low frequencies to — ► 
the second term in the brackets in Eq. (12) diverges in any dimension 
d ^ 2. This means that in a disordered film the quantum correction to 




(12) 
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the classical conductivity diverges and this signals (but, of course, does 
not prove) the possibility of the localization. 

The scaling theory of the localization and the discovery of the new 
diffusion modes was revolutionary in the theory of disordered systems 
but how can these findings be related to the Wigner-Dyson theory? 

Actually, Eq. (12) is the key to understanding that there can be 
something beyond the classical Eq. (10) in small metal particles. One 
needs only realizing that the case of small metal grains corresponds to 
the zero dimensionality of the integral in Eq. (12). Due to the finite size 
the values of the momentum k are quantized such that k a = 2irn a /L a , 

where n a = 0, ±1, ±2, ±3 , a = (x,y,z), and L a is the size of the 

grain in the a -direction. 

At low frequencies u> <C Dq/L 2 the most important contribution in 
the integral in Eq. (12) comes from the zero harmonics with k = 
and one can see that this contribution strongly diverges when u> — > 0. 
Moreover, the quantum correction is proportional to A/uj, where A is 
given by Eq. (11) and this is what one can expect from Eqs. (4-6). 

So, Eq. (12) really signals that something nontrivial can happen in 
metal grains and the WD theory is not excluded. The diffusion modes 
play a prominent role and it seems, at first glance, that one should 
merely write proper diagrams and sum their contribution. However, 
even if this were possible this would hardly correspond to Eq. (4-6). 
The problem is that the expansion in terms of A/u cannot take into 
account the oscillating part in Eqs. (4-6) even in principle. Summing the 
diagrams one can hope to reproduce only non-oscillating asymptotics of 
these equations. This means that another approach has to be developed. 

The possibility to demonstrate that the energy level and wave func- 
tions statistics can really be described by the RMT came first with the 
development of the supersymmetry approach [13, 14]. This method is 
based on a representation of Green functions of a disordered metal in 
terms of an integral over both commuting and anticommuting variables. 
Singling out excitations with the lowest energy (diffusion modes) one 
can reduce calculations to a supermatrix non-linear model. The super- 
symmetry method allowed to prove for the first time that the level-level 
correlation function for disordered particles is really described by Eqs. 
(4-6). Later, using the supersymmetry technique Verbaarschot, Wei- 
denmiiller and Zirnbauer [15] have derived Eqs. (4-6) starting directly 
from the Gaussian ensembles, Eq. (1). 

When deriving the non-linear u-model for metallic particles it was 
very important that they contained disorder. However, it is not the nec- 
essary condition. Several years later Bohigas, Gianonni and Schmidt [16] 
conjectured that the RMT should describe correctly spectral properties 
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of quantum systems which are chaotic in their classical limit. In par- 
ticular, the Wigner-Dyson statistics had to be observed in clean metal 
particles (quantum billiards) provided their shape was such that classi- 
cal motion would be chaotic. Their hypothesis was made on the basis 
of extensive numerics. For a review of the subsequent activity in these 
fields the book [17] is a good reference. 

Historically, the description of disordered systems with a non-linear 
(T-model has been suggested by Wegner [18] using the replica method 
and integration over conventional complex variables. In the first work 
there were problems with convergence of functional integrals and there- 
fore the replica approach was further developed in the publications [19] 
and [20] . The a- model of Ref . [20] was obtained by the integration over 
conventional variables and, as a result, the group of the matrices Q was 
non-compact. In contrast, the starting point of Ref. [19] was a repre- 
sentation of Green functions in terms of integrals over anticommuting 
(Grassmann) variables and this lead to a compact group of the matrices 
Q. 

Although both the replica and supersymmetry approach are equiva- 
lent when doing the perturbation theory in the diffusion modes, the lat- 
ter method is much more efficient for non-perturbative calculations like 
the study of the level-level correlations. This had become clear shortly 
after the works [19, 20] were finished and this drawback of the replica ap- 
proach motivated the development of the supersymmetry one. It should 
be noticed that recently the oscillating behavior of the level-level correla- 
tion function R (uj) has been obtained [21-23] using the compact replica 
cr-model of Ref. [19]. However, the procedures used in these references 
are considerably more complicated than the calculations by the super- 
symmetry method and the limit of low frequencies u < A is still hardly 
achievable. 

It is fair to say that the replica approach allows including electron- 
electron interactions in a comparatively easy way [24, 25] and this has 
been the main motivation in the attempts [21-23] to obtain non-perturbative 
results within the replica technique. At the same time, it was believed 
for a long time that an inclusion of the electron-electron interaction into 
the supersymmetry scheme was impossible. However, this is not quite 
so and, at least, not very strong interaction can be incorporated in the 
supermatrix a- model [26]. 

It follows from this discussion that the supersymmetry approach is 
better suitable for making connections with the RMT and therefore the 
present lectures contain discussions based on this method only. It is pos- 
sible neither review here all works made in this direction nor present all 
details of the calculations. For a more detailed information see the book 
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[4] and more recent reviews [27-31]. It is relevant to mention here that 
the word supersymmetry has appeared in the condensed matter physics 
in the publication by Parisi and Sourlas [32], who discovered a complex 
symmetry in a model describing ferromagnets in a random magnetic 
field. They used a concept of superspace including both commuting and 
anticommuting variables. 

Two other related directions of the use of the RMT are reviewed at 
this school by Boris Altshuler (Quantum Chaos) and Jac Verbaarschot 



In the next sections I want to present the main ideas of the super- 
symmetry approach and show how it can be used for both the level cor- 
relations and wave functions statistics. It will be shown how to obtain 
the Wigner-Dyson statistics and how to go beyond it. A new devel- 
opment concerning a generalization of the supersymmetric <T-model to 
more complicated situations will be outlined in the last section. 

1. Supersymmetry method. 
Super mat hematics 

The supersymmetry method is based on the use of the so called Grass- 
mann variables Xii * = 1> 2,..., n (the elements of the Grassmann alge- 
bra) that are introduced in a completely formal way. These are abstract 
objects but in many cases abstract mathematical constructions drasti- 
cally influence the development of physics. For example, nobody can 
dispute the usefulness of complex numbers for physics, but what is the 
physical meaning of sj— 1? Here I want to remind the reader basic 
formulae concerning definitions and operations with objects containing 
combinations of the Grassmann variables and conventional numbers (su- 
permathematics) . 

The Grassmann variables are some mathematical objects obeying the 
following anticommutation rules [33] 



for any 1 < i, j < n. 

The anticommutation rules, Eq. (13) hold in particular for i = j and 
we see that the square of an arbitrary variable \i ls zero 



For any anticommuting variable x one can introduce its "complex con- 
jugate) x* ■ It is assumed by the definition that (x*)* = ~X, sucn that 
the "square of the modulus is "real" 



(QCD). 



{Xi,Xj} = XiXj + XjXi = 



(13) 




(14) 



\XiXi) XiXi XiXi 



(15) 
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The anticommuting variables Xii Eq. (13-15), remained not very use- 
ful until Berezin introduced integrals over these variables. The integrals 
are nothing more than formal symbols introduced as follows 

J ' d Xl = j d X * = 0, J XidXi = J XidXi = 1 (16) 

It is implied that the "differentials" d%i, d X * anticommute with each 
other and with the variables Xii Xi 

Wxi,dXj} = {dXi,dx*j} = {d X *,d X j} = (17) 

{dXi,Xj} = {dXi,Xj} = {dx*,Xj} = {dx*,Xj} = 

The definition, Eq. (16), is sufficient for introducing integrals of an 
arbitrary function. If such a functions depends only on one variable Xi it 
must be linear in X i because already X \ = 0- Assuming that the integral 
of a sum of two functions equals the sum of the integrals we calculate the 
integral of the sum with Eq. (16). The repeated integrals are implied by 
integrals over several variables. This enables us to calculate the integral 
of a function of an arbitrary number of variables. 

The most important for the development of the supersymmetry method 
are Gaussian integrals. The direct integration shows that the following 
relation is fulfilled 



„ n 

J exp {- X + A X ) H d X *d Xi = DetA (18) 



where A is an n x n Hermitian matrix and 
/ Xi \ 



X 



X2 



\Xn ) 



X + ={X\ Xl ■ ■ Xn) (19) 



Eq. (18) differs from the corresponding equation for the commuting 
variables by giving det A instead (det Ay 1 . This remarkable difference 
is the basis of the supersymmetry method presented in these lectures. 
In addition to Eq. (18) one can write one more useful integral 

T JXiX* k exp(- X +A X )ll? = idxtd Xl , ^ 

Je W (- X +A X )U? =1 dxld Xl [Ji >* [ZU) 

In contrast to the integral /, Eq. (18), the integral I2, Eq. (20), is com- 
pletely similar to the corresponding integral over conventional numbers. 
Eq. (20) can be proven by the differentiation of In I in A^. 
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The next step is the introduction of supervectors and supermatrices. 
An n + m component supervector is introduced as 



X 

s 



where the n- component vector \ is defined in Eq. (19). 
component vector S has a similar form 

S 2 



\s m ) 



(21) 
The m- 



(22) 



but its components are conventional complex numbers. 

In analogy with conventional vectors one can introduce the Hermitian 
conjugation 

$+ = (<D T )* (23) 

and the scalar product 



(24) 



a=l 



a=l 



A linear transformation F in the space of the supervectors converts a 
supervector $ into another supervector $ 



(25) 



F = 



(26) 



Of course, the supervector $ must have the same structure, Eq. (21), 
as the supervector <J>. This imposes a restriction on the structure of the 
supermatrix F corresponding to the linear transformation F : it has to 
be of the form 

a a 
p b 

In Eq. (26) a and b are n x n and m x m matrices containing only 
commuting variables, a and p are n x m and m x n matrices consisting 
of anticommuting ones. Matrices having the structure, Eq. (26) can be 
called supermatrices. 

Two supermatrices F and G of the rank (m + n) x (n + m) are as- 
sumed to multiply according to the conventional rules 



m+n 



(27) 



z=i 
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and one can see that FG is a supermatrix of the same form. In order 
to define the supertranspose F T of the supermatrix F one should use 
the notion of the scalar product of two supervectors, Eq. (24). Again, 
by analogy with the conventional definition the supermatrix F T is in- 
troduced as 

$fF T $ 2 = (F$i) T $ 2 (28) 

The transpose of a conventional matrix is obtained by transposing its 
indexes. This is not as simple for the supermatrices. Writing out the 
scalar product on both sides of Eq. (28) explicitly and using the anti- 
commutation relation, Eq. (13), one can see that the supermatrix F T is 
equal to 

F-=(fr /) (29) 

where a T , b T ', a T , and p T stand for the conventional transposition of the 
matrices a, b, a, and p. 

Using the scalar product, Eq. (28) one obtains immediately 

(iqF 2 ) T = if if (30) 

The Hermitian conjugate F + of the matrix F can be defined in a stan- 
dard way 

F+ = (F T )* (31) 
Combining Eqs. (15, 30, 31) one can obtain 

(F 1 F 2 ) + = if if and (F+) + = F (32) 

The latter equality shows that the operation of the Hermitian conjuga- 
tion is inverse to itself. The same is not generally true for the transpo- 
sition 

(F T f + F (33) 

A very important operation in the theory of conventional matrices is 
taking the trace of a matrix. If one takes the trace of a product of sev- 
eral matrices it is invariant under cyclic permutations of the matrices. 
However, due to the presence of anticommuting elements a proper op- 
eration for the supermatrices should be defined in a different way. The 
supertrace STrF of matrix of the form, Eq.(26) is defined as 

STrF = Tra - Trb (34) 

where the symbol Tr stands for the conventional trace. 

Although somewhat strange, the definition, Eq. (34) is very useful 
because it is this operation that provides the invariance under the cyclic 
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permutations. We obtain for arbitrary supermatrices Fj of the form, Eq. 
(26) 

STrFiFi = STrF 2 F 1 (35) 

and 

STr (FiF 2 ...F„) = STr {F^F^.F^) (36) 

In addition to the supertrace it is convenient to introduce a superdeter- 
minant of the supermatrix F 

In SDetF = STr In F (37) 

The super determinant SDetF can also be written as 

SDetF = Det (a - ab' 1 p) Detb' 1 (38) 

The connection between the superdeterminant and supertrace enables 
us to prove immediately the multipliticity of the superdeterminant 

SDet (FiF 2 ) = {SDetFx) (SDetF 2 ) (39) 

The rules of the operations with the supervectors and supermatrices are 
very convenient because they are similar to those of conventional linear 
algebra. In fact, one can manipulate superobjects in exactly the same 
way as conventional objects. This simplifies calculations with quantities 
containing both types of variables considerably. 

Now we can write Gaussian integrals over supervectors that generalize 
the integrals over conventional complex numbers or Grassmann variable. 
A direct calculation shows that 

I s = I exp d$*d$> = SDetF, (40) 



= n d X *dxi n ds*dSi 

i=l k=l 

and 

= /$^exp(-$+F$)d$*d<S = /p-i\ (41) 

2 / exp d<S>*d$ [ )ik K ' 

The formulae written in this subsection give complete information 
about integrals over the Grassmann variables, supervectors, and super- 
matrices. This information will be directly used for constructing the 
supersymmetry method. 
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Physical quantities as integrals over supervectors. 
Averaging over disorder. 

Eq. (41) enables us to express physical quantities in terms of func- 
tional integrals over supervectors. The form of the integrals that will 
be obtained is such that averaging over disorder can be performed the 
beginning of all calculations. 

I start with the Schrodinger equation for electrons without any electron- 
electron interactions but in a presence of an external potential containing 
both regular and irregular parts. The regular part can describe potential 
walls and other features of the system whereas the irregular part H\ of 
the Hamiltonian stands for disorder. The Schrodinger equation takes 
the form 

H4> k = e k <f> k , H = H + H 1 , (Fx) = (42) 

where 4> k and e k are eigenf unctions and eigenvalues, respectively. The 
angular brackets (...) stand for the averaging over disorder. 

The most important physical quantities can be expressed in terms 
of retarded G^ and advanced G A Green functions of the Schrodinger 
equation. Using the spectral expansion the Green functions G^' A can 
be written in the form 

G n,A (r; r>) = ^ ^» = ^ G n^ k (r) ^ {r>) (43) 

k k 

These functions satisfy the equation 

(e - H) G?' A (r, r') = S (r - r') (44) 

The average density of states p (e) (this quantity is proportional to 
A -1 , Eq. (2)) can be written as 

(p (s,r)) = {^2 <t>l ( r ) fa (r) S (e - e k )^ (45) 
= -(lmG^(r,r)) (46) 

7T 

whereas the level-level correlation function R (u) takes the form 



R(u) = 




We see from Eqs. (45, 47) that, as soon as we are able to average 
the Green functions or their products over the disorder, the quantities 
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of interest are found. However, this cannot be done directly using Eqs. 
(45, 47) and we need another representation for the Green functions. 
Of course, one can do perturbation theory in the disorder potential but, 
as it has been already discussed, such an approach can hardly help in 
obtaining the Wigner-Dyson statistics. 

What will be done now is writing the Green functions in a form that 
would be suitable for averaging over the disorder in the very beginning. 
This can be conveniently done with the integrals over the supervectors. I 
want to present here the main scheme only. All necessary details can be 
found in the book [4]. As the main interest is to calculate the level- level 
correlation function R (uj) , Eq. (47), all formulae will be written for this 
case. Calculating the density of states, Eq. (45), is a simpler and less 
interesting task. 

Let us introduce 8-component supervectors tp consisting of 4-component 
supervectors t/' 1 and ip 2 such that 

lb m _( m \ 1 ( X m *\ rn_ 1 ( S m * \ 

Y> -{ v m),0 ~ 73 ^ x m ),V - -j= y sm J, (48) 

m = 1,2. 

For the supervectors 4> of the form of Eq. (48) one can define, in 
addition to transposition and Hermitian conjugation, the operation of 
the "charge conjugation" 

ijj = [Ctp) T , 4> m = ( e m , v m ) (49) 

In Eq. (49), T stands for transposition, and C is the supermatrix of the 
form 

QVtin j^mn I ] 

\ c 2 J 

where A is the diagonal supermatrix 




with 1 the unity 4x4 unity matrix. 
The matrices c\ and C2 have the form 

Cl = ( 1 ~0 ) ' C2 = ( 1 I ) (51) 

The function R (uj) is determined by products G A G R , G R G R and G A G A . 
The last two products are not interesting because the average of their 
products is equal for weak disorder to the product of their average. When 
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calculating these quantities the diffusion modes discussed in the Intro- 
duction do not appear. Therefore let us concentrate on the calculation 
of the product G A G R . Using Eqs. 40, 41) we can write this quantity as 

Gt. (r, r) Gf (r', r') = (^) / ^ ( r ) P ( r ) ^ ¥ ( r « X P (-L) Zty 

- " ' (52) 
where the Lagrangian L has the form 

L = i J i> (r) ^-# - C/ (r) - ^ (w + z<5) V (r) dr (53) 

where U (r) the impurity potential. The operator Hq equals 

H = e (-iV r ) - £ + | (54) 

where e (— iV r ) is the spectrum. 

I would like to draw attention at this point that the weight denomina- 
tor is absent in Eq. (52), which contrasts the analogous integrals, Eqs. 
(20,41). Actually, this is a consequence of the fact that F in Eq. (41) is 
taken as unity in the space of the 2x2 supermatrices when writing Eq. 
(52). In other words, the weight denominator is absent due to the differ- 
ence of the results of the Gaussian integration over the anticommuting 
variables %, Eq. (18), (one obtains DetA) and the corresponding formula 
for the integration over conventional numbers that gives (DetA)^ 1 . 

The possibility to write the Green function in the form of Eq. (52) 
without the weight denominator is the reason why the integration over 
the supervectors is used. As the weight denominator is absent one can 
immediately average over impurities in Eq. (52). 

Let us assume for simplicity that the distribution of the random po- 
tential U (r) in Eq. (53) is specified by Eq. (8). Then, the averaging 
over the random potential is simple and one obtains again Eq. (52) but 
now the Lagrangian L should be written as 



L 



dr (55) 



Eq. (55) shows that we have reduced the initial disordered problem to 
a regular model with an if; 4 interaction. Of course, it does not help to 
solve it exactly but now we can use approximations well developed in 
theory of interacting particles. 

The Lagrangian L in Eq. (55) is similar to those studied in field theory. 
Let us remark that at uj = the Lagrangian is invariant under rotations 
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of the supervectors in the superspace because it depends on the square 
of "length" only. The frequency u violates this symmetry and, if one 
uses an analogy with spin models, plays the role of an "external field". 
The violation of supersymmetry by the frequency is due to the fact that 
Eq. (55) is written for the product G R G A (the presence of the matrix 
A is a direct consequence of this) . The symmetry would not be violated 
in the corresponding integrals for (G R G R ^ and (G A G A ^ . The averaging 
of the simpler Lagrangian corresponding to the density of states (one 
averaged Green function) results also in a model with the interaction 
ip 4 , but the supersymmetry in this case is not violated. The diffusion 
modes discussed previously exist only as a result of the violation of the 
supersymmetry (Goldstone modes). 

Spontaneous breaking of the symmetry and 
Goldstone modes. Non-linear cr-model. 

It is clearly not possible to calculate any correlation function with 
the Lagrangian L, Eq. (55), exactly (except for the case of the one- 
dimensional chain or the Bethe lattice, where one can write recurrence 
equations. Further calculations will be performed in the limit of large 
mean free times r, which correspond to a weak interaction in the La- 
grangian L. However, the use of the standard perturbation theory as 
we have seen is impossible even in this limit because of the existence 
of the diffusion modes and so one should try to use non-perturbative 
approaches. 

One of the standard approaches used for such a type of theories is the 
mean field approximation. According to this scheme one simplifies the 
V> 4 interaction replacing pairs tptp by their averages. For the interaction 
V> 4 there can be six different pairings, which can be written as follows 




(56) 
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In Eqs.(56), the symbol 

<->e// = / (...)exp(-L e// )^ 

stands for the functional averaging with the effective Lagrangian L e f f = 
Lq + L\ + L2 + L3, where Lq is the quadratic part of the Lagrangian L 
in Eq. (55), and a and (5 stand for the components of the supervectors 
tp and if). 

In fact the terms in L2 and L3 are equal to each other. The average 
(V'aV'a) renormalizes the energy e and is not important. The averages 
in L2 and L3 can be both commuting and anticommuting variables, 
depending on the subscripts a and (3. The final Lagrangian L e ff takes 
the form 



tpdr 



(57) 



with the 8x8 supermatrix Q satisfying the following self-consistency 
equation 

Q = h We// (58) 
Calculating the Gaussian integral in Eq. (58) with the help of Eq.(41) 
we obtain 

Q = v„l^$f < 59 > 

90 (p) = 1 [e ( P ) + -e + -(u + i6)A+^] (60) 

The integral over the momenta p in Eq. (60) has both a real and an 
imaginary part. As concerns the imaginary part the main contribution 
comes from the region \e (p) —e\ 3> r _1 , lo and, therefore, is proportional 
to the unit matrix 1. This contribution leads to a small renormalization 
of the energy e. Assuming that this energy has already been renormalized 
we can forget about the imaginary part of Q and concentrate on the real 
part. The main contribution to the real part of Eq. (60) comes from the 
region \e (p) — e\ ~ r _1 <^ e. Introducing the variable ^ = e (p) — e we 
can rewrite Eq. (60) in the form 

Q= uiH (ij+<5)A+ sr* <6i> 

Eq. (61) determines the contribution to the real part of Q only and has 
at lo / one evident solution 

Q = A, lo + (62) 
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However, putting ui = we see that any Q of the form 

Q = VAV (63) 

where V is an arbitrary unitary supermatrix, VV = 1, satisfies Eq. (61). 
Here the symbol of conjugation " — " for an arbitrary matrix A means 
the following 

A = CA T C T (64) 

with C from Eqs. (49, 51). The supermatrix Q is self conjugate, Q = Q. 

The degeneracy of the ground state, Eq. (63), leads to the existence of 
the low-lying Goldstone modes, and their contribution to physical quan- 
tities should be taken into account properly. These Goldstone modes are 
just the diffusion modes discussed in the Introduction. In the language 
of spin models the degeneracy of the solution, Eq. (63) is equivalent to 
the degeneracy due to an arbitrary spin direction. 

Substituting the mean field solution, Eq. (63), into the effective La- 
grangian into Eq. (57) we obtain zero, which shows that this approxima- 
tion is not sufficient. In order to describe a contribution of the diffusion 
modes we have to take into account slow variations of the supermatrix 
Q in space. This can be done assuming that the supermatrix V is a slow 
function of the coordinates r and expanding in the gradients of V (or, 
equivalently, Q). 

As a result of the expansion in the gradients and small frequencies 
one can obtain a functional F [Q] describing slow variations of the su- 
permatrix Q in space 

f [Q] = ^ STr J [ D o ( V( 2) 2 + 2i ( w + i5 ) K ®\ dv ( 65 ) 

where -Do is the classical diffusion coefficient, Eq. (7), and the superma- 
trix Q is described by Eq. (63). The free energy functional, Eq. (65), 
has the form of a non-linear u-model. 

In order to calculate, e.g. the level-level correlation function R (u>) , 
Eq. (47) , one should express it in terms of a functional integral over the 
supermatrix Q. As a result, this function takes the form 

R ( UJ ) = l-^F2 Re j Qu W Qn ( r ex P (" F IQ]) DQdvdv' (66) 

where V is the volume of the system. In Eq. (66) the superscripts of 
Q enumerate "retarded-advanced blocks" , the subscripts relate to the 
matrix elements within these blocks. 

Eqs. (65, 66) show that the calculation of physical quantities for 
disordered systems can be reduced to study of a supermatrix non-linear 
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cr-model. There is a variety of physical problems that can be solved by 
considering this model in different dimensionalities. To some extent, the 
problem of the level statistics in a limited volume is the simplest one 
because it corresponds to the zero-dimensional cr-model (no dependence 
of Q on the coordinates). 

Adding magnetic or spin orbit interactions results in additional "ex- 
ternal fields" partially breaking the symmetry. If these interactions are 
not very weak they simply cut some degrees of freedom. As a result, one 
comes again to Eq. (65) but with a reduced symmetry of the supermatrix 



There can be three different types of the symmetries: 

1. In the absence of both magnetic and spin-orbital interactions the 
system is invariant under time reversal and spatial inversion. It will be 
called "orthogonal" because, as we will see, one comes in small particles 
to the Wigner-Dyson statistic for the orthogonal ensemble. 

2. In the presence of magnetic interactions the time reversal symmetry 
is violated and this will be called unitary ensemble. 

3. In the absence of magnetic interactions but in the presence of 
spin-orbital ones, one obtains the symplectic ensemble. 

These three ensembles lead to quite different results in different situ- 
ations and each of them should be considered separately. 

Level statistics in a limited volume. 

Eq. (63) specifies a general form for the 8x8 supermatrix Q. In other 
words, it obeys the constraint 



However, this constraint is not sufficient to determine unambiguously 
the precise structure of Q because Eq. (67) could correspond both to ro- 
tations on a sphere and on a hyperboloid. It turns out that the symmetry 
of the supermatrix Q is more complex than those corresponding to these 
two possibilities. The supermatrix Q consists of two parts: one describ- 
ing rotations on the sphere and the other on the hyperboloid, such that 
the group of rotations of Q is a mixture of compact and non-compact 
groups of rotations [4] . These two parts are glued by anticommuting el- 
ements. We can describe the explicit form of the supermatrix Q writing 
it in the form 



Q- 



Q 2 = l 



(67) 



Q 



UQ U, 



(68) 



where 




(69) 
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and uu = 1, vv = 1. All Grassmann variables are included in the 4x4 
supermatrices u and v. These matrices contain also some phases. Their 
form is simple and not very important for our discussion. More details 
can be found in Ref. [4]. 

The form of Qq is more interesting and can be written as 

Qo= {-ism§ -cose )' e -{ o e 22 ) (70) 

The 2x2 matrices are different for different classes of the symmetry. 

We see that the symmetry of Qo corresponds to a group rotations on 
both the sphere and the hyperboloid. The explicit form of the matrices 
0n and 622 can be expressed as follows 





M l ! q g | • "22 



n "2 
h 0i 



with < 9 < 7T, 9i > 0, 02 > for the orthogonal ensemble, 

\ * Y 0i 



11 — n a i "22 



o / ' V o 0i 

with < < 7T, 0i > for the unitary ensemble, and 



(71) 



(72) 



n "2 

'11-1 a a I i "22 



6 ° ' (73) 



02 01 / ' V 

for the symplectic one. 

Eqs. (65,66, 70-73) specify the non-linear supermatrix er-model. What 
remains to do for the level-level correlation function, Eq. (66), is to cal- 
culate the functional integral over Q. Many of the physical quantities can 
also be expressed in a form of a correlation function of the supermatrices 
Q with the free energy functional F [Q] , Eq. (65). 

Let us consider now the level-level correlation function in a limited 
volume. The functional integral, Eqs. (65, 66) can be considerably 
simplified in the limit of small frequencies. In a finite volume one can 
expand the supermatrix Q in Fourier series. Then, it is not difficult to 
understand from Eq. (65) that for uj <C Dq/L 2 , where L is the sample 
size, only the zero space harmonics is essential. In this case one can 
integrate over the supermatrices Q not varying in space and the integral 
for the function R(u>), Eq. (66), becomes just a definite integral over 
several variables. This integral takes the following form 

R( L o) = ^-^Rej QllQll exp (-F [Q]) dQ (74) 
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where Fq [Q] takes the form 

FolQ] = 'JLl^lsTr(A Q) (75) 

One can say that Eqs. (74, 75) determine a zero dimensional non-linear 
<7-model. In general, it is clear that the dimensionality of the cr-model 
is determined at not very high temperatures by the geometry of the 
sample. For example, the one dimensional a-model describes a long 
wire of a finite thickness. 

In order to calculate the integral over Q in Eq. (74) it is very con- 
venient to use Eqs. (70, 73). We see immediately that the supermatrix 
U, Eq. (69), drops out from Fq [Q] entering the pre-exponential in Eq. 
(74) only. This allows us to integrate first over the elements of U, and 
reduce the integral to the variables On, # 22 only. The integration over 
the Grassmann variables is, according to Eq. (16), a very simple task. 
Actually, one has to calculate also Jacobians arising when changing from 
the integration over Q to the integration over the "eigenvalues" A, Ai, A2. 
As a result, one comes to the following integrals for all three ensembles 

f°° f°° f 1 (AiA 2 -A) 2 (l-A 2 ) , N 

R orth (co = 1 + Re / / / V 76 

Ji Ji J-i (A 2 + A 2 + A 2 -2AAiA 2 -l) 

x exp [i (x + iS) (A1A2 — A)] d\id\2(l\ 



] f'OO f'l 

Runit (oj) = 1 + - Re J J exp [i (x + iS) (Ai - A)] dAdAi (77) 

f°° f 1 f 1 (A-AiA 2 ) 2 (A 2 -l) , N 

Ji Jo J-i (A 2 + Af + A^ - 2AiA 2 A - 1) 

x exp [i (x + id) (A — A1A2)] dAi(iA 2 dA 

We see that the integrals over the supermatrix Q can be reduced to 
integrals over the "eigenvalues" A, Ai, A 2 . Depending on the ensemble 
one obtains twofold or threefold integrals. Calculation of such integrals 
is a much simpler task than calculation of integrals over a large num- 
ber TV of variables encountered in RMT[2]. Although the reduction to 
the integrals over the eigenvalues has been carried out for the level-level 
correlation function only, the corresponding manipulations for studing 
other physical quantities are the same. The only thing that remains 
to be done when calculating different physical quantities is to write a 
proper pre-exponential and carry out integration over the elements of 
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supermatrices u and v entering the pre-exponential only. Then one ob- 
tains integrals over the variables A, Ai, A2 analogous to those in Eqs. 
(76, 77, 78). 

The integration over A and Ai in Eq. (77) for the unitary ensemble is 
simple. At first glance, the integrals for the orthogonal and symplectic 
ensembles look terrifying. However, they can be calculated by Fourier 
transforming the integrals from the frequencies u to the real time t. As 
a result, one comes to the Wigner-Dyson formulae, Eqs. (4-6), which 
demonstrates that the level-level correlation function for a small metal 
particle is really the same as that given by the RMT. This is how one 
proves the relevance of the RMT for disordered systems [13]. Actually, 
the model of the disordered metal was the first microscopic model for 
which the Wigner-Dyson statistics had been proven. 

Of course, this is possible for not very high frequencies u> <C Dq/L 2 . 
In the opposite limit uj S> Dq/L 2 , the situation is no longer zero di- 
mensional but one can do perturbation theory in diffusion modes. This 
calculation was done in Ref. [34]. 

One can also come to the zero dimensional cr-model starting from the 
Gaussian distribution for the random matrices, Eq. (1), in the limit of 
large N . This was done in the review, Ref. [15]. Therefore, the super- 
symmetry method can be considered as an alternative to the method of 
the orthogonal polynomials [2] in RMT. 

Using the non- linear cr-model, Eq. (65), one can consider thick dis- 
ordered wires. This corresponds to the one dimensional cr-model. In 
this case one can use a transfer matrix technique that allows to reduce 
calculation of a one dimensional functional integral to solving of an ef- 
fective "Schrodinger equation" . The exact solution found for this model 
[35] proves localization of all states (vanishing of the conductivity) for 
any arbitrarily weak disorder. In the language of random matrices this 
quasi-one-dimensional model corresponds to a model of random banded 
matrices [36]. 

The exact solution can also be found for a model with disordered 
grains connected in a such a way that they constitute a Bethe lattice. 
For this model, using recursion relations one can write a non-linear in- 
tegral equation[37]. It was demonstrated that within the model on the 
Bethe lattice one could have a metal-insulator transition with a very un- 
usual critical behavior. The model on the Bethe lattice has been shown 
to be equivalent to models of certain sparse random matrices [38] . Sparse 
matrices are relevant for description of diluted spin models, some com- 
binatorial optimization problems [39] and other interesting systems. 

Properties of two dimensional disordered metals can be studied using 
a renormalization group scheme [18-20, 13, 14] 
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2. Wave functions fluctuations in a finite 
volume. Multifractality. 

General definitions. 

In this Section statistics of wave functions is discussed. Investigation 
of wave functions is complimentary to study of energy levels. In the 
language of random matrices, energy levels correspond to eigenvalues of 
the matrices whereas the wave functions relate to eigenvectors. 

Study of wave functions has become popular in condensed matter 
physics not long ago with the development of mesoscopic physics. One 
can study, e.g., electron tunneling through so called quantum dots, which 
are actually small quantum wells. At low temperatures one can reach a 
resonance tunnelling regime when the electron tunnel via one resonance 
level. In this case the tunnelling amplitude is very sensitive to the wave 
function of the resonance state. 

Wave functions can also be measured in optical and acoustic res- 
onators where they are electromagnetic or sound waves, respectively. 

We start with the standard Schrodinger equation 

H<p a (r) = e a 4> a (r) (79) 

that determines the eigenenergies e a and eigenfunctions 4> a (r) . 

We assume that a finite volume V is considered, such that the spec- 
trum of the energies e a is discrete. 

The full statistics of the wave functions (p a (r) at a given point r can 
be described by the following distribution function / 

f (*) = A (e 6 (* - mi 2 ) 6 ( £ - ( 8 °) 

The function / (t) , Eq. (80) , gives the probability that the square of the 
absolute value of the wave function (intensity) at the point r and energy 
e is equal to t. The distribution function / (t) and the wave functions 
(j) a (r) are assumed to be properly normalized such that 

to = 1, h = V- 1 (81) 

where t n (V) are coefficients of the so called inverse participation ratio 




These coefficients indicate very sensitively the degree of localization of 
states through their dependence t n (V) on the volume of the system. In 
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a pure metal or a ballistic chaotic box where the wave functions extend 
over the whole system one has 

t n oc V' n (83) 

If disorder makes the localization length L c , at which the typical wave 
functions decay, is much shorter than the sample size L ~ V l / d , the co- 
efficients t n are insensitive to L. However, a very interesting information 
about the development of localization can be gained through an analysis 
of t n (V) for small samples with L < L c . 

As soon as the localization length L c exceeds the sample size, any 
length scale disappears and, in the language of the coefficients t n , this 
is described as 

Vt n oc L- T{n \ t (n) = (n - 1) d* (n) (84) 

where d* (n) may differ from the physical dimension d of the system 
and be a function of n. This function gives the values of the fractal 
dimensions d* (n) for each n. 

If the behavior of the wave function is described by Eq. (83), the 
fractal dimension d* coincides with the physical dimension d. We will 
see that the fractal dimension of a system obeying the RMT coincides 
with the physical dimension. In such a situation, although the amplitude 
fluctuations are possible, they are not very strong. 

Once we assume that the envelope of a typical wave function at a 
length scale shorter than L c obeys a power law 4> (r) oc with a single 
fixed exponent /x < d/2, the set of the coefficients t n reveals d* = d — 2\i 
different from d but the same for all n > dj (2/x) . This is when one 
speaks of fractal behavior with the fractal dimension d*. 

If d* (n) is not a constant, that signals a more sophisticated structure 
of the wave functions. They can be imagined as splashes of multiply 
interfering waves at different scales and with various amplitudes, and 
possibly, self-similarity characterized by a relation between the ampli- 
tude t of the local splash of the wave function and the exponent n (t) of 
the envelope of its extended power law tail. 

We will see below that the multifractality of the wave functions of 
two dimensional weakly disordered conductors is the most general prop- 
erty of these systems as soon as the sample size L does not exceed the 
localization length L c [40]. 

Porter-Thomas distribution. 

Before starting more complicated calculations I would like to present 
here what is known about the distribution of wave functions from nuclear 
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physics (see, e.g. Refs. [41, 42]) where it was studied for description 
of level width fluctuations in neutron scattering. The wave function 
fluctuations are obtained there again from the RMT. 

To start the calculation one should choose and arbitrary basis of eigen- 
functions p m (r) and expand the function 4> n (r) in this basis 



oo 



<P n (r) = V- 1 ' 2 a n m p m (r) (85) 

m=l 

where V is the volume. It is convenient to truncate the basis to a finite 
iV-dimensional set and take the limit N — ► oo, as is usually done in 
the RMT. The main statistical hypothesis is that all coefficients ci m n 
are uniformly distributed. The only restriction on the coefficients {a mn } 
is imposed by normalization of the wave functions, and the probability 
density P ({a mn }) can be written as 



P({a mn }) = ^(Y,\^n\ 2 ~l\ (86) 

\m=l / 



where Q n is the solid angle in N dimensions. Because of the truncation 
of the basis the condition of completeness of the basis {p m } should be 
written in the form 

oo 

J> 2 „(r)^| 2 = iV (87) 

m=l 

where ff is an iV-dimensional vector with components {p m } . The dis- 
tribution function of the intensities W (v) at the point r is introduced 

as 

W (v) = ^- J S(v- \aff\ 2 ^j 5 (|a| 2 - l) da (88) 

where the vector a is an iV-dimensional vector with components {a mn } . 
The distribution function W (v) , Eq. (88) is defined in such a way that 



/ 



W (v) dv = 1 (89) 



In the unitary ensemble, one should integrate over complex vectors a. 
Integrating first over the component of the vector a parallel to the vector 
fj and using Eq. (87) one obtains 

"V-mzJ'iM-i 1 -*))" (90> 
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where a± is the component perpendicular to ff. The remaining integra- 
tion in Eq. (90) can be carried out easily. Taking the limit N — > oo one 
obtains for the unitary ensemble a simple formula 



Computing the integral in Eq. (88) for real vectors a and ff one can 
obtain the distribution function for the orthogonal ensemble 



The functions W (v), Eqs. (91,92) satisfy the normalization conditions, 
Eq. (89). The distribution functions W (v) are universal and do not 
depend on details of models for disorder. The amplitudes v are related 
to t from the previous section as v = Vt. The functions W and / are 
related to each other accordingly. Eqs. (91, 92) are usually referred to 
as the Porter-Thomas distribution. 

Non-linear cr-model and the statistics of wave 
functions. 

Now we concentrate on the calculation of the distribution function 
/ (t) , Eq. (80). At first glance, this task does not seem easy. In the pre- 
vious section we were able to reduce the level-level correlation function 
to a functional integral over 8x8 supermatrices Q. It became possi- 
ble because the level- level correlation function R(lo) , Eq. (3), could be 
expressed in terms of the product of two Green functions, Eq. (47). 
Actually, the size 8 x 8 of the supermatrices is determined by the fact 
that only two Green functions are needed. 

So, in order to calculate the distribution function / (t) , Eq. (80), we 
have to make two necessary steps: 

1. To express Eq. (80) in terms of the Green functions. 

2. To express products of the Green functions in terms of an integral 
over supervectors ip. If we really want to make explicit calculations the 
supervectors ip should not have too many components. 

It turns out that both the steps are possible and the necessary number 
of the components of the supervector ip is just 8 for the orthogonal 
ensemble and 4 for the unitary one. 

The step 1 is done introducing Green functions G^ A for a system 
with smeared levels 



W (v) = exp (—v) 



(91) 




(92) 




(93) 
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Then, Eq. (80) can be written as 

/ («) = A ^ lim £ S (t - ^G* (r, r)) 5 ( e - £a ) ^ (94) 

= ^^&(/ 5 (*-T G ^ (r ' r) ) (^(r',r')-G^(r',r')))dr' 

In Eq. (94) one should first take the limit /3 — > and then 7^0. Since 
the distribution function / (t) is represented in terms of a function of only 
two Green functions at two points r and r/ one can express it in terms of 
an integral over 8-component supervectors ip using the Wick theorem. 
The corresponding Lagrangian is the same as the one for the level-level 
correlation function, Eq. (55), provided one replaces the frequency uj by 
the level width 7. 

The derivation of the corresponding <r-model is standard and one 
comes to the following expression for the distribution function / (t) 

/(i) = lim| j ' STr[i^Q{T))5[t-^STr[i^Q{T ))) (95) 
dr 

X (-F[Q})DQ— 
where the free energy functional F [Q] has the form 



TTU f 

~8J 



STr 



Do (VQ) 2 - 7AQ (r) (96) 



and r D is the "observation point" . When the system can be described by 
the zero dimensional a- model the distribution function / (t) does not 

depend on t . However, beyond the 0D approximation, this function can 

(1 2) 

also be a function of the coordinates. The matrices ir b in Eq. (95) 
select from the supermatrix Q its boson-boson sector and have the form 

(i) ( vr 6 \ (2) ( \ /0 0\ 

< ={ J' < ={o n)' n= {o 1 J (97) 

As we have discussed previously, the a- model is noncompact, Eqs. (70- 
73). Therefore, in order to avoid divergent integrals over the variables 
#22; we must calculate the integrals keeping 7 finite. The limit 7^0 
can be taken only at the end of the calculations. However, it is not 
very convenient to keep an additional free parameter, and it is better 
to get rid of the parameter 7 at an earlier stage. This can be done by 
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integrating over the zero space harmonics of Q in the very beginning of 
the calculations. 

To carry out this procedure one should represent the supermatrix 
Q (r) in the form of Eq. (63) and change the variables of integration 
V (r) to V (r) as V (r) = V (r ) V (r) . This leads to supermatrices Q : 
Q (r) = V (r„) Q (r) V (r ) . In terms of the new variables V (r) and Q (r) 
the gradient term in Eq. (96) preserves its form, but now the condition 

V(v ) = 1, Q(r G ) =A (98) 

has to be fulfilled. Changing the variables of the integration for all points 
r ^ r G from Q (r) to Q (r) one obtains a new free energy functional that 
does not contain V (r G ) or Q (r G ) . These variables enter only the pre- 
exponential and the term with 7, and, hence, the integral over V (r D ) 
can be computed without making approximations. The result of the in- 
tegration contains only the variables Q (r) with the boundary condition, 
Eq. (98). This means that the reduced a- model obtained in this way 
operates only with relative variations of the field Q with respect to its 
value at the observation point. 

The limit 7^0 simplifies the computation because the main contri- 
bution to the integral over the variable B\ entering the parametrization, 
Eqs. (70-73), is from cosh#i D ~ I/7 (for simplicity we are considering 
the unitary ensemble but the final results are similar for all ensembles). 
After standard manipulations one can express the distribution function 
f (t) in the form 

^) = ^^§^' / exp(-F[Q,*l)l>g(r) (99) 

V dt JQ(r )=A V L J/ 

where the free energy F[Q,t] has the following form 



F[Q,t] = \ J STr itvDq (vq) 2 - 2t AIIQ 



dr (100) 



The matrix II selects from Q its noncompact "boson-boson" sector. 

If t is not very large one can take into account the zero space harmonics 
of Q only. Taking into account Eq. (98) we can just put everywhere 
Q = A, which leads us immediately to the Porter-Thomas distribution, 
Eq. (91) (and Eq. (92) for the orthogonal ensemble). 

Nonetheless, this would be only an approximate procedure because the 
value Q (r) = A does not correspond to the minimum of the functional 
F when t 7^ 0. An equation for the minimum can be found by taking into 
account the noncompact variable 6\ under the condition at the boundary 
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and at the observation point 

nV0i = 0, 6»i (r ) = (101) 

where n is the unit vector at the boundary and perpendicular to it. 

Using Eq. (101) one writes the equation for the extremum solution 9 t 
in the form 

A r t = --^-ex P (-0 t (r)) (102) 

where A r is the Laplacian. The solution 6 t (r) of Eq. (102) has to be 
substituted into the energy functional F, Eq. (100), which takes the 
form 

F t = \j [™D {V6 t ) 2 + 2t exp (-6 t )] dv (103) 

It is remarkable that Eq. (102) for the non-trivial vacuum of the reduced 
u-model in two dimensions is exactly the Liouville equation known in 
the conformal theory of 2D quantum gravity [43, 44]. Within this model 
one has to calculate the functional integral over all 6 with the free energy 
functional determined by Eq. (103). Although this can lead to helpful 
analogies [45], results that might be anticipated in this way can be used 
only as intermediate asymptotics. 

Eq. (102) for the minimum looks very similar to a saddle point equa- 
tion derived in Ref. [46] when considering the problem of long-living cur- 
rent relaxation. However, the non-linear term is different, which leads 
to different solutions. 

The most interesting is the solution of Eq. (102) in two dimensions 
where it can be found exactly. However, the exact solution is somewhat 
cumbersome and I write here its asymptotics at distances r much smaller 
than the sample size L (but exceeding the mean free path /) 

exp (-d t ) « (l/rf» (104) 

where \x is a parameter depending on disorder. 

With the same accuracy, the free energy of the vacuum state can be 
approximated by 

F t « 4tt 2 uD {fi + n 2 In (L/l)} (105) 
The parameter fj, can be determined from the following equation 

a « where ze =T = = y ' ' (106 

p 21n(L/i) 2ir 2 vD V ' 

In principle, Eqs. (104-106) determine the distribution function / (t) 
for arbitrary t (p, must remain small, though). However, analytical ex- 
pressions can be written only in the limiting cases T<1 and T » 1. In 
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these limits, the distribution function / (t) , Eq. (80), can be written as 
[40] 

f exp(-Vt[l- | + ...]), T<1 

exp(-^^T), T»l W 

where A is a normalization constant. 

We see from Eq. (107) that at small values of the amplitudes, such 
that T <C 1, the distribution function f (t) agrees with the Porter- 
Thomas distribution, Eq. (91), thus proving the latter for disordered 
systems. In this limit one can make expansion in T[47]. 

At large f (T > 1) the function / (t) has log-normal asymptotics that 
is strikingly similar to the asymptotics of the distribution function of the 
local density of states or conductances discovered by Altshuler, Kravtsov 
and Lerner [48] who came to this result considering renormalization of 
terms high gradients in the cr-model. Even the numerical coefficients 
in the exponentials are the same, although, of course, the logarithms 
contain different variables. It appears that the log-normal form is really 
universal. The slower decay of the distribution function / (t) at large t 
is due to localization effect. Unfortunately, until now it is not clear how 
the growth of the high gradient terms is related to the existence of the 
non-trivial vacuum considered in this Section. 

As concerns the coefficients of the inverse participation ratio t n , Eq. 
(82), they show the multifractal behavior, Eq. (84). Using Eqs. (102,106) 
one comes to the following expression for the fractal dimension d* (n) 

d*(n) = 2-n (A^vDoy 1 (108) 

We see that even for a weak disorder the fractal dimension d* (n) strongly 
deviates from 2. Of course, it cannot become negative and this deter- 
mines the region of the applicability of Eq. (108). 

3. Recent and possible future developments. 

In the preceding sections it was demonstrated how one can come to 
the random matrix theory starting from a model of a disordered metallic 
particle. This became possible using the supersymmetry method. Actu- 
ally, to the best of my knowledge, the model of a disordered metal was 
the first microscopic model for which the Wigner-Dyson statistics was 
confirmed. 

Starting from the first works [7, 13] where the relevance of the Wigner- 
Dyson theory to the disordered systems was suggested and proven, a 
huge number of problems have been attacked using these ideas. Cal- 
culations were performed either assuming that RMT was applicable for 
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the description of small particles (they are often called "quantum dots" ) 
and using methods of the RMT [2] or making direct computation start- 
ing from a disordered metal and applying the non-linear supermatrix 
<7-model. Reviewing all these application within several lectures is im- 
possible even though several related topics are considered in this volume 
by Boris Altshuler and Jac Verbaarschot. At this point I can only refer 
again to recent reviews [4, 6, 17, 27-31, 49] and apologize in case if some 
references are missing here. 

The selection of topics of the present lectures was motivated mainly by 
the desire to give a feeling of how to calculate within the supersymmetry 
method both level and wave function correlations. We have seen that 
one could obtain results that agreed in a certain region of parameters 
with the predictions of the RMT and, at the same time, go beyond the 
RMT. 

Essential conditions for the derivation of the cr-model were the absence 
of the electron-electron interaction and a sufficiently high concentration 
of impurities. For the problem of the level statistics, the latter condition 
corresponds to the case when the mean free path I is much smaller than 
the sample size. In other words, an electron can scatter many times on 
the impurities in the bulk before it reaches a boundary of the sample. 

At the same time, the RMT was initially suggested for description of 
complex nuclei, where disorder is absent but interactions are strong. A 
natural question that can be asked is: Can one prove the relevance of 
the RMT for clean or/and interacting systems? Clearly, the supermatrix 
cx-model discussed in the previous sections is not applicable in these 
situations and one needs a generalization of this method. 

It seems that really new ideas are necessary in order to achieve this 
goal. Nevertheless, first steps towards constructing more general schemes 
have been done and I want to present here the main ideas of the new 
approaches. 

Supersymmetry with interaction. 

From the beginning of the use of the supersymmetry method it was 
clear that the method could be applicable for non-interacting particles 
only. The method is based on the result of the Gaussian integration, Eq. 
(18), that gives DetA instead of the usual {DetA)~ l . Introducing an in- 
teraction results in non-quadratic terms in the Lagrangian. Therefore 
the trick with writing Green functions in terms of a Gaussian integral 
without a weight denominator does not work anymore. This is the rea- 
son why, in contrast to the replica approach where a proper cr-model 
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has been derived long ago [24, 25], introducing an interaction into the 
supersymmetry scheme was believed to be impossible. 

To some extent, it is true and it is not clear how to include the in- 
teraction into the supersymmetry exactly. However, a weak interaction 
can really be included without considerable difficulties [26]. The initial 
electron model with the interaction is not super symmetric and cannot 
be made super symmetric by a transformation. This is why one cannot 
get rid off the weight denominator. 

The main idea of Ref. [26] is to replace approximately the initial 
electron model by an effective super symmetric model. This is possible 
for any disorder in the limit of a weak interaction. The effective model 
takes into account the most important (Hartree-Fock type) diagrams for 
any fixed configuration of impurities. Then, the derivation of the proper 
cx-model is quite standard. 

The resulting non-linear cr-model resembles very much the replica 
model of Finkelstein[24, 25] but contains supermatrices and does not 
have replica indices. The supermatrices Q contain, in addition to those 
for the non- interacting systems, indices for Matsubara frequencies. One 
should also write properly spin indices. As a result, the supermatrix 
non-linear cr-model for electron systems with interaction takes the form 
(unitary ensemble) 



where 71 and 72 are certain operators acting on the supermatrices Q, 
and Ti and T2 are scattering amplitudes characterizing the interaction 
(they are different from those obtained for the replica cr-model [24, 25]. 
As usual, one has the constraint Q 2 = 1 but now the product of two 
supermatrices includes summation also summation over the Matsubara 
frequencies. 

Using the sigma-model Eq. (109) renormalization group equations of 
Refs. ([24, 25]) have been reproduced in the first order in the interaction 
constants and there is a hope to use this model for non-perturbative 
calculations. 

Method of quasiclassical Green functions. 

Although the interaction is included in Eq. (109), it is written for 
systems with considerably high concentration of short range impurities. 
Studying problems for clean systems or systems with a long range disor- 
der one needs a different scheme. Statistical properties of clean chaotic 
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systems are covered in this school by Boris Altshuler and I do not review 
them here. Instead, I want to concentrate on calculational schemes. 

The saddle point equation (61) is not a good approximation for clean 
systems and systems with long range disorder and we cannot follow the 
same way as the one used for diffusive models. The method that I want 
to present now is based on using quasiclassical Green functions. 

Introducing an 8 x 8 matrix function G (r, r') as 

G(r,r') = 2^(r)^(r')) v , (110) 
we can write in a standard way the following equation for this function 



H 0r + U(r)+A^A + iJ(r) 



G (r, r') = id (r - r') (111) 



where the subscript r of Hq t means that the operator acts on r. The 
notations are the same as in Eqs. (53, 54), and J(r) is a source term 
that allows to extract more complicated correlation functions. 

Conjugating Eq.(lll)we obtain another equation for the matrix G(r, r') 
with the operator H 0r > acting on its second variable 



G(r;r') 



H 0r , + U(r')+A^±^ + iJ(r' 



iS(r - r') (112) 



Until now no approximations have been done and Eqs. (Ill, 112) are 
exact. Now we can use the assumption that the potential U (r) changes 
slowly on the wavelength A^. If the mean free path I for the scattering on 
the random potential exceeds A f the Green function varies as a function 
of r — r' at distances of the order of \f but, at the same time, is a 
slow function of (r + r') /2. The Fourier transform G p ((r + r') /2) of G 
(r, r') respective to r - r' has a sharp maximum near the Fermi surface. 
In order to cancel large terms we subtract Eq. (112) from Eq. (111). 
Using the assumption that the potential U (r) is smooth and expanding 
it in gradients we obtain in the lowest order 



G p (R)+^^[A,G p (R)] +i [J(R),G p (R)] = 

(113) 

where R = (r + r') /2 and [, ] stands for the commutator. When deriving 
Eq. (113), not only the potential U (r) but also the function J (r) was 
assumed to be smooth. 

The dependence of the Green function G p (R) on |p| is more sharp 
than on other variables. In order to avoid this sharp dependence we 
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integrate Eq. (113) over |p|. Of course, this procedure makes a sense 
for very large samples when the level discreteness can be neglected. 

The most interesting contribution in the integral over |p| comes from 
the vicinity of the Fermi-surface. A contribution given by momenta 
considerably different from pp is proportional to the unity matrix an 
drops out from Eq. (113). 

Introducing the function g n (r) 

g n (r) = lj G pn (r) d£, £ = (114) 

where n is a unite vector pointing a direction on the Fermi surface, we 
obtain the final quasiclassical equation 

{v F nV - V f V r U (r) dn) <?„ (r) + ' + [A, 9n ( r )] _ [J, g n ] = 

(115) 

where 

d 

d n = V n - n, V n = -[n x nx — 

an 

Eq. (115) should be complemented by a boundary condition at the 
surface of the sample. Considering a closed sample we assume that the 
current across the border is equal to zero. This leads the boundary 
condition at the surface 

9n± ( r )l surface = 9-a ± ( r )\ surface (H6) 

where is the component of the vector n perpendicular to the surface. 

Eq. (115) is similar to an Eilenberger equation written long ago in 
superconductivity theory [50]. As in the theory of superconductivity, 
the solution for the Eq.(115) satisfies the condition 5„(r) = 1. Eq. (115) 
is written for a non-averaged potential U (r) and it is valid also in the 
absence of the long range potential. 

The quasiclassical equation, Eq. (115), has been written first by 
Muzykantskii and Khmelnitskii, Ref. [51] who guessed a functional for 
which Eq. (115) is just a condition for an extremum. Then, they pro- 
ceeded to work with this functional without estimating fluctuations near 
this minimum. 

It came later as a surprise that, actually, one could write the proper 
solution of Eq. (115) in terms of a functional integral over supermatrices 
exactly [52]. The exact solution for Eq. (115) can be written as 
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g n (r) = Z-\J\ [ Q n (r) exp (-^j[Q n (r)]) DQ n , 
$[Q n (r)] = Str J drdn[Af n (r)(v F nV r 

-p^V r E/(r)V n )T n (r) + ^ + ^ A - J(r)) Q n (r)], (117) 
Qn(r) = T n (r)Af n (r), T„(r)T n (r) = 1 
In Eq. (117), the partition function Z[J (r)] is 

Z[J}= [ exp (-^[Q n (r)}) DQ n (118) 

and the integration is performed over the self-conjugate supermatrices 
Qn = Qn (r) satisfying the following relation 

QoW = i (H9) 

We see that the quasiclassical Green function g n (r) can be written in 
the form of a functional integral over supermatrices Q n (r) depending 
both on the coordinates r and the direction of the momentum n and 
satisfying the constraint, Eq. (119). The first term in the free energy 
functional is written in terms of the supermatrices T n rather than Q n . . 
However, it can be written in a form of a Wess-Zumino-Novikov-Witten 
term containing the supermatrices Q only [51]. Writing this term one 
should introduce an additional coordinate varying at the interval [0, 1] . 

The model described by the functional <£, Eqs. (117), is usually re- 
ferred to as a "ballistic a- model" . The partition function Z [J] , Eq. 
(118), is unity at J = due to the supersymmetry and this allows us to 
average (if necessary) over the smooth potential U (r). 

The method of quasiclassical Green functions suggested here for a 
static external potential does not seem to be restricted by the non- 
interacting case. There are indications that it can be generalized to 
describe clean interacting systems. Of course, in this case one should 
write the quasiclassical equations in time representation because the in- 
teraction mixes states with different energies. Study of interacting sys- 
tems with this method may be a very interesting direction of research. 

As concerns attempts to prove the Wigner-Dyson statistics for clean 
non- interacting systems one can try to start with the functional Eq. 
(117). At first glance, we should simply restrict ourselves with the in- 
tegration over Q depending neither on the coordinates, nor on the mo- 
menta. Then, the functional $ would contain only the last term and we 
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would have the zero dimensional cr-model, which leads immediately to 
the WD statistics. 

However, a very important question is whether one averages over the 
potential U (r) or puts U (r) = and averages over the spectrum. In 
the former case one gets after averaging over U (r) an additional term 
in functional $ quadratic in gradients. This term leads eventually to a 
suppression of non-zero harmonics and one can really obtain the zero- 
dimensional a- model (see [52] and references therein). 

The situation with U (r) = and averaging over the energy is more 
interesting. Everything depends on whether the system is classically 
integrable or chaotic. It is just the situation for which the authors of 
Ref. [16] made their hypothesis. 

It turns out that within the model with the functional Eq. (117), 
and U (r) = one cannot come to the zero-dimensional cr-model. There 
is a common consensus that a "regularizer" (see e.g. [53]) containing 
something like square of gradients in coordinates or momenta is neces- 
sary in the correct ballistic cr-model. Aleiner and Larkin [54] argued that 
in order to come to the zero dimensional cr-model one had to take into 
account diffraction, which is clearly absent in the ballistic cr-model, Eq. 
(117). They did not manage to include the diffraction in their calcula- 
tional scheme microscopically and mimiced it by introducing artificial 
quantum impurities that would correspond to the potential U (r) in Eq. 
(117). This allowed them to come to the zero dimensional cr-model, 
confirm the WD statistics and calculate corrections to it. It is worth 
emphasizing that the effective potential U (r) was very weak such that 
the computation was done in the ballistic regime. 

As concerns the real physical diffraction, it cannot be directly included 
in the cr-model using the quasiclassical scheme and a more sophisticated 
approach is necessary. 

Beyond the quasiclassics 

We have seen that the solution of the equations for the Green func- 
tions, Eqs. (Ill, 112), can be written in the quasiclassical approximation 
in terms of the functional integral over 8x8 supermatrices Q n (r) . For 
certain problems this approximation is not sufficient and the natural 
question is: can we do better than that and find a solution for the Green 
functions in terms of a functional integral over supermatrices valid at all 
distances including those of the order of the electron wavelength Xp ? 

This attempt has been undertaken recently in Ref. [55], where an 
integral of such a type was suggested for a solution of Eq. (111). The 
idea is rather close to the one known in field theory where it is called 
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bosonization[56]. The final expressions obtained in Ref. ([56]) are rather 
complicated and this method, to the best of my knowledge, has not 
evolved into an efficient calculational tool. 

However, the supersymmetric form of the Green functions considered 
here seems to promise more and the derivation is rather simple. I follow 
here a simpler derivation of Ref. [57]. 

What I want to show now is that the Green function, Eq. (Ill), 
can be represented exactly as an integral over supermatrices Q (r, r') 
depending on two coordinates r and r/ 



G{v,v') = Z- 1 [J] j Q(r,r')exp(-<I>[Q])DQ 
where Z[J] is a new partition function 

Z[J] = J exp(-$[Q]) J DQ 
and the functional $>[Q] has the form 



(120) 



(121) 



<D[Q] = -Str / [H 0r + U(r) + 



to + iS 



A 



+ 



x 5(r - r')Q(r,r')drdr' 
istrlnQ-^Str J J(r, r')Q(r', r)drdr' 



(122) 



where J (r, r') is a source term. The structure of the supermatrix Q(r, r') 
in the integral Eq.(120) should be the same as that of the Green function 
G(r, r'), i.e. be the same as of the product ip(r)ip(r'). In particular, this 
means that Q(r, r') is self-conjugated 



Q(r,r / )^CQ T (r',r)C T = Q(r,r') 
In order to prove Eq.(120) we write the following identity 

8 exp (— iStrlnQ) 



(123) 



- hz-\j\ 



ex P I ~2 Str 



<SQ(r",r 

to + i5 



A + iJ 



-Q(v",r')dr" 



Q DQ 



iS(r - r') 



(124) 



and integrate over Q by parts. The derivative 5/5Q should act now on 
both Q and the exponential. At this point, the supersymmetry plays 
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a crucial role. Differentiating first the supermatrix Q we obtain the 
supermatrix product (5/5Q)Q. As the number of the anticommuting 
variables in the sum over the matrix elements is equal to the number of 
the boson ones and the derivatives have the opposite signs, this matrix 
product vanishes. Differentiating the exponential only we come to the 
following equation 

Z- 1 [J] J dr" (A 0r + U (r) A + ijj (r, r") x 

J Q(r",r')exp(-<S>[Q])DQ = i5(r-r f ) (125) 

Eq. (125) proves immediately that the integral Eq.(120) does satisfy 
Eq.(lll) and we have really the alternative representation of the Green 
function in terms of an integral over the supermatrices Q. 

Making the Fourier transform Q (r', r") in the difference r'— r" (Wigner 
transformation) one can express the functional [Q] in terms of the 
variables Q p (r) , where p is the momentum and r is the center of mass 
(r' + r") /2. Then the free energy functional <I>, Eq. (122) can be written 
as 

®[Q} = \ Jstr [Hj(x) *Q(x) -ilnQ(x)]dx (126) 
where x = (r, p) is the coordinate in the phase space, 

Wj(x) = H (p) + U (r) + + iJ (p, r) (127) 

is classical Hamilton function. 

The product * of two matrices A (x) and B (x) is defined by Moyal 
formula 

4 I V r Vp-V r Vp ) 

A(x)*B(x) = A(x)e 2 \ VB(x). 

The scheme of calculations using the Wigner representations and the 
star product "*" is known as Weyl symbol calculus [58]. This method is 
convenient for quasiclassical expansions. 

If the potential U (r) is smooth, one can simplify Eq. (126) and come 
again to Eq. (117). This procedure is described in Ref. [55]. The func- 
tional <I>, Eq. (126) has a form of the Lagrangian of a non-commutative 
field theory [58]. The method suggested here can naturally be called 
"superbosonization" . 

At the end I have to warn that Eqs. (120-122, 126) are not complete 
yet because nothing has been said about the contour of integration over 
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the supermatrices Q. This was not very important in Refs. ([55, 57]) 
where a saddle point approximation was used. However, generally, this 
question requires a more careful study and this is a subject of a current 
work. In case if the difficulties are overcome, the superbosonization can 
become very useful in the theory of random matrices for more general 
models than the one described by the conventional Eq. (1). 

4. Summary 

In these lectures I tried to achieve two goals: 

1. The random matrix theory is to a large extent a phenomenological 
theory. Therefore, it is very important to have examples when it can be 
obtained starting from a microscopic model. The model of the disordered 
metals considered here is the first one for which the relevance of the 
RMT has been proven. The proof became possible with the help of the 
supermatrix non-linear a-model first derived for other purposes. 

2. Having presented the derivation of the a-model I demonstrated 
how one obtains the Wigner-Dyson statistics from its zero-dimensional 
version. However, in many situations the supersymmetry method allows 
to go beyond the Wigner-Dyson model and obtain completely different 
results like the log-normal distribution of the amplitudes of the wave 
functions of Section 2. Of course, as soon as one has a Hamiltonian 
one comes to random matrices. However, generally it is not clear how 
to write the distribution function for these matrices in each particular 
situation and the supersymmetry method can play an important role for 
investigation of microscopic models. 

I wanted by no means to oppose the RMT and the supersymmetry 
method to each other. They can be considered as complimentary meth- 
ods, although with a considerable overlap. We see at this school that 
the random matrices find more and more applications in many fields of 
physics, which is an exciting development. I believe, in many cases the 
supersymmetry technique can also be useful in these new applications 
and one should keep in mind a possibility of using this scheme. 
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